Stretching of polymers in a random three-dimensional flow 
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Behavior of a dilute polymer solution in a random three-dimensional flow with an average shear 
is studied experimentally. Polymer contribution to the shear stress is found to be more than two 
orders of magnitude higher than in a laminar shear flow. The results indicate that the polymer 
molecules get strongly stretched by the random motion of the fluid. 

PACS numbers: 83.10 Nn, 83.50 Ws 



Behavior of flexible polymer molecules in a solution 
in various flows is one of the basic problems of polymer 
physics 0. Flows that have been studied most exten- 
sively are simple shear and unilateral extension. Those 
flows are rather easy to create in the laboratory and the 
experimental data are most straightforward to interpret. 
As a result of many mechanical [QJ|] , optical Jl|,^,0 , theo- 
retical d|| and especially the most recent single molecule 
investigations 0^), the major issues of polymer dynam- 
ics in these two types of flow seem to be resolved now. In 
contrast to it, there are virtually no experimental data 
on dynamics of polymer molecules in a generic case of a 
complex flow with a random velocity field. In this let- 
ter we present results of our measurements of mechanical 
properties of a dilute polymer solution in a random three- 
dimensional (3D) flow. These results directly indicate 
that polymer molecules can become strongly stretched 
in such a flow. 

Stretching of polymer molecules manifests macroscopi- 
cally as growth of optical anisotropy of the polymer so- 
lution and growth of the stress, r p , which is due to poly- 
mers [0. By definition, t p = r — t s , where r is the total 
stress in the solution, and t s is the stress in the same 
flow of a pure solvent. 

Three-dimensional flows with a random velocity field is 
quite a general class of flows, which also includes turbu- 
lence. Such flows usually occur at high Reynolds number, 
Re, that implies large fluid velocity and big size of the 
tank. The high Re flows are usually non-homogeneous 
in space even in statistical sense. Therefore, bulk opti- 
cal measurements are difficult to interpret, and so are 
measurements of mechanical stresses, a, at a wall of the 
flow system. In a random, high Re flow the momentum 
is transferred across the fluid not only by the molecu- 
lar forces contributing to the stress, t, but also by the 
macroscopic motion of the fluid. The latter produces the 
Reynolds stress p < ViVj >, where p is the fluid density, 
and v is fluctuating part of the velocity V. Thus, we have 
a momentum balance equation for the average values of 
the stresses 



+ T P + P < ViVj >, 



(1) 



where the right hand side is taken in the bulk of the fluid. 
That is, in order to calculate the average of r p , one needs 
to evaluate the average of p < ViVj > over the whole vol- 
ume of the flow, which is rather difficult to carry out. 
A way to surmount the problem of the high Reynolds 
stress is to create a random flow with a low Reynolds 
number. Such a situation is realized in the case of the 
elastic turbulence, which we have reported recently ||. 
It is an apparently turbulent flow that can arise in poly- 
mer solutions as a result of non-linear elastic effects at 
arbitrarily low Re. 

We carried out our experiments in a swirling flow between 
two parallel plates. Polymer solution was held in a sta- 
tionary cylindrical cup with a flat bottom (lower plate). 
A coaxial rotating upper plate was just touching the sur- 
face of the fluid. The cup was mounted on top of a com- 
mercial rheometer (AR-1000 of TA-instruments) with a 
normal force transducer. The upper plate was attached 
to the shaft of the rheometer, which allowed precise con- 
trol of its rotation velocity, f2, and measurements of the 
torque, T, and the average stress, a, applied to it. The 
sidewalls of the cup were made transparent which allowed 
measurements of the flow velocity in the horizontal plane 
by a laser Doppler velocimeter, LDV. In the first exper- 
iment, the radii of the upper plate and of the cup were 
R = 30 mm and i?2 = 43.6 mm, respectively, and the 
distance between the plates was d — 30 mm. Because of 
the big gap and short distance to the side wall, the shear 
rate profile was rather non-homogeneous even in a lam- 
inar flow. So, we defined the average applied shear rate 
as the ratio between a in a laminar flow, and the fluid 
viscosity. It was j a v — 2.78£lR/d. The whole set-up was 
put into a transparent box, and temperature of fluid was 
stabilized to better than 0.05 °C by throughflow of air. 
We used 25 ppm solutions of high molecular weight 
polyacrylamide (M,„=18,000,000, Polysciences) in vis- 
cous sugar syrups. The syrups were made out of sucrose 
and sorbitol (corn sugar) in a proportion of 1 : 2, and 1% 
of NaCl was added to stabilize the ionic contents. In the 
first experiment the total sugar concentration was 76.3%, 
and the solvent viscosity at the temperature of the ex- 
periment, 18 °C, was 77 s = 1.36 Pa-s. The polymer part 
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of the solution viscosity, rf = (t — T s )/j, was chang- 
ing at j av , that we applied, in the range from 0.12?7 S 
to 0.08?7 S , decreasing with the shear rate. At such low 
values of rf /rj s effect of overlapping of separate polymer 
molecules is supposed to be minor. 

When flexible polymers get stretched in a shear flow, they 
also get aligned along the flow direction Q . This leads to 
a difference in the normal stresses in the streamwise and 
transverse directions, that for a flow between two rotat- 
ing plates is N\ = — t zz . Here (r, </>, z) are cylindrical 
coordinates. This first normal stress difference results 
in a normal force, F n , which pushes the plates apart 
(the "hoop stress") jij. So, N± can be evaluated from 
measurements of F n . For flexible polymer molecules, 
Ni is connected with the polymer relaxation time, A, by 
Ni = 2rf\ A f 2 [Q. For our solution we measured A = 6.3 
s. 

In the first experiment we evaluated dependence of the 
stress at the upper plate, cr, and the normal force, F n , on 
the Weissenberg number, Wi — Xj a v (The role of the 
Weissenberg number in the elastic turbulence is similar to 
the role of the Reynolds number in the usual turbulence 
in normal fluids ||.) Simultaneously, one velocity com- 
ponent was measured in the center of the set-up. Figure 
la, curve 1, shows a divided by the stress, <Ti am , mea- 
sured in a laminar flow with the same applied shear rate, 
as a function of Wi. It resembles very much Fig. 2 from 
Ref.10. The flow at low Wi is laminar. At Wi of about 5 
transition from the laminar flow to the elastic turbulence 
occurs. It manifests in growth of the flow resistance and 
in onset of fluctuations of the fluid velocity, Fig. lb. The 
amplitude of the velocity fluctuations increases with Wi. 
The ratio a/ai am increases with Wi as well, reaching a 
value of about 13 at Wi = 22. The Reynolds number, 
Re = flRdp/r] s , is only 1.3 at the highest Wi. 
As we showed elsewhere the flow of a polymer solu- 
tion at high Wi (above the transition) bears all features 
of developed turbulence. The flow velocity changes ran- 
domly in space and in time, and the fluid motion is ex- 
cited in a broad range of spatial and temporal scales. The 
subject of the current study is stretching of the polymer 
molecules in this random 3D flow. The elastic turbulence 
itself is driven by the polymer stresses, which are gener- 
ated by the stretched polymer molecules. However, ori- 
gin of forces driving a flow is not directly relevant to the 
problem of polymer stretching. Extension of polymers 
can only depend on local properties of the flow velocity 
field along the trajectory of the fluid element, which con- 
tains the polymer molecules. 

In spite of turbulent character of the fluid motion at high 
Wi, the Reynolds stresses were so small, that the corre- 
sponding term in Eq.(l) could be totally neglected. So, 
at Wi=22 characteristic amplitude of the velocity fluc- 
tuations, V rms , measured in different points was about 
0.7 mm/s (see Fig. lb for V rms in the center). Thus, the 
Reynolds stress could be estimated as pV 2 ms = 7 • 10~ 4 



Pa, while the stress a was 91 Pa. Eq. (1) is written in 
a general form valid for momentum transfer in the direc- 
tion perpendicular to a wall. In our case we have cir- 
cular symmetry. So, it is torque, T, at the upper plate, 
which is measured, and it is angular momentum, which 
is transferred from the plate through the liquid to the 
stationary cup. The average flux of angular momentum 
due to the solvent stress, t 3 , is defined by r) s r 2 (doj/dr) 
and r} s r{dCo/d4>), where u> = V^/r is local average angu- 
lar velocity of the fluid. The turbulent flow changes the 
distribution of to compared to the laminar case, leading 
to larger gradients of ui near the upper plate [|l0| . Nev- 
ertheless, the boundary conditions on w at the surfaces 
of the upper plate and of the cup remain the same. So, 
for our estimates we take the volume average of the flux 
of the angular momentum due to solvent to be the same 
as in the laminar flow. Then the whole increment in T 
and a should be solely due to growth of average t p in the 
fluid bulk. Taking contribution of the solvent stresses to 
a to be the same as in the laminar flow, we get for the 
contribution of the polymer stresses to a the curve 2 in 
Fig. la. This curve gives the ratio of the average shear 
stresses due to polymers in the turbulent and the laminar 
flows, T p /r^ m . This ratio reaches a value of 170, which 
is an evidence of strong stretching of the polymers in the 
turbulent flow. 

Dependencies of average N\ on Wi measured in the elas- 
tic turbulent flow and in a laminar flow are shown in 
Fig.lc. One can learn that Ni in the turbulent flow be- 
comes about an order of magnitude higher than in the 
laminar flow with the same j a v This is another evidence 
of stretching of the polymers by the random 3D fluid mo- 
tion. 

A specific feature of the elastic stresses is that they do not 
turn to zero immediately after the fluid motion stops, but 
rather decay with their characteristic relaxation time, A 
0. That is how they can be clearly distinguished from 
the viscous stresses, which decay instantaneously. So, a 
way to independently examine the origin of stresses ap- 
plied to the upper plate is to stop its rotation and to 
measure decay of the stresses. This was the objective 
of our second experiment, Fig. 2. It required a higher 
polymer relaxation time. Therefore, the concentration 
of sugars was increased to 80.3%, that gave 77 s of about 
7.2 Pa-s and A of about 30 sec at the temperature of the 
experiment, 15 °C. The size of the set-up was reduced by 
a factor of 2, R = 15 mm, R2 = 21.8 mm, d — 15 mm. 
In this small set-up the dependence of a on Wi for the 
polymer solution from the first experiment was the same 
as in Fig. la (see also Ref.10). The characteristic torques 
were 8 times smaller, however, that enabled a sharp stop 
of the upper plate. 

Rotation of the upper plate was started abruptly at an- 
gular velocity ft = 0.07 s _1 corresponding to Wi ~ 6 and 
Re pa 3 • 10 -3 . The rotation continued for a while, and 
then, at a time moment taken as zero, ft was abruptly 
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brought to zero. In the first run, curve 1 in Fig. 2b, the 
time of the rotation was short, about 115 seconds. It was 
just enough for the polymers to get properly stretched by 
the primary shear and for the polymer stress to reach its 
laminar flow value, t? . So, by the time the rotation 
was stopped at t = 0, the major part of torque was due 
to the solvent, T s w 10 /iN-m, and it relaxed almost im- 
mediately. (Characteristic viscous diffusion time is very 
small, d 2 p/r] s — 0.04 s.) This torque value, 10 /iN-m, 
was also measured just after rotation had been started 
and polymers had not had time to get stretched yet. Ex- 
trapolation of the curve 1 in Fig. 2b to t = +0 gives the 
value of the torque due to polymers in the laminar flow, 
Tf am = 0.85 M N-m. 

In the second run the time of rotation, about 4500 sec, 
was long enough for the transition to developed elastic 
turbulence to complete. The torque reached its satu- 
rated value, T = 82.5 /iN-m, corresponding to applied 
Wi. From Fig. 2a one can see that the fluid velocity was 
fluctuating during the time, that the torque increased. 
As a result of the turbulent fluid motion, the slowly re- 
laxing part of the torque increased by almost two orders 
of magnitude. One can learn from the inset that the 
relaxing torque can be reasonably extrapolated to 72.5 
/iN-m at t — +0. Thus, the immediately relaxing part 
of the torque is again due to the solvent shear stresses, 
T s « 10 /iN-m, and the results of these relaxational mea- 
surements are quite consistent with the suggestion, that 
the whole increment in the torque is due to growth of 
the polymer contribution. If inertial effects were impor- 
tant, they would obviously lead to torques of the oppo- 
site sign. The polymer torque T p = 72.5 /iN-m at t = 
means average increase of t p by 85 times compared with 
the laminar flow. 

Our measurements show that as a result of a secondary 
random 3D flow, superimposed on the primary applied 
shear, the polymer contribution to shear stress, r^, can 
become as much as 170 times larger. If linear elasticity 
of the polymer molecules is assumed, the stress is related 
to the polymer extension via jf, ~< R z R,f, >, where 

R is the vector connecting edges of a polymer chain |Q . 
Then the factor of additional extension of the polymer 
molecules due to the random 3D flow is about 13. In 
fact, measurements of the stress in the turbulent regime 
can be viewed as a mechanical test, where shear is applied 
to a polymer solution, which flows turbulently with zero 
average shear rate. Then, one can say that polymer con- 
tribution to viscosity, defined as rf = T^/jav, increases 
by more that two orders of magnitude because of the 
turbulent flow. In theories of polymer dynamics rf is al- 
ways connected with an effective volume occupied by the 
polymer chains [El. Thus, our experimental results sug- 
gest that the effective volume can significantly increase, 
and polymer molecules can get strongly stretched in a 
random 3D flow with zero averages. 



Dynamics of a polymer molecule in a random 3D flow 
were first considered by Lumley p[ , and have been re- 
vised recently jl2|,[l3) ■ It is suggested that the flow is al- 
ways homogeneous on the scale of a polymer molecule, so 
that velocity field in some vicinity of a molecule with the 
center at r is given by V — V(r ) + n ■ (r — r ). Possible 
complex structure properties of the flow at larger scales 
are not important for the issue of the polymer stretching. 
It only depends on statistics of the tensor of rate of de- 
formation, k, which varies randomly in time and space. 
If the flow is truly 3D, k always has an eigenvalue with a 
positive real part, so that there is a direction along which 
pure extension occurs [§El]. The direction and the rate 
of extension change randomly as a fluid element is rotat- 
ing in the flow and moving along its trajectory. Neverthe- 
less, if k remains correlated within finite time intervals, 
the overall statistical result of such random motion will 
be exponential divergence of two closely spaced material 
points. In a turbulent flow an estimate for the correlation 
time of k is given by the inverse of the velocity gradients 
themselves. Then the average Lyapunov exponent, a, for 
divergence of two material points is given by the rms of 
longitudinal velocity gradients with a prefactor of order 
unity. Thus, in the statistical sense, a random 3D flow 
acts as an extensional flow with dV x /dx — a, where the 
direction, x, of the maximal extension of a fluid element 
is changing randomly in time and space Jll],[l2| ■ In such a 
flow polymer molecules should become vastly stretched, 
if the condition aX > 1/2 is fulfilled §,0,EJ,EJ|. In fact, 
strong polymer stretching in a turbulent flow is expected 
to occur even at smaller average a Il2j . 
A probability distribution function (PDF) of a longitu- 
dinal velocity gradient, dV$/(rd4>), for Wi = 12 in the 
big set-up (first experiment) is shown in Fig. 3. The rms 
of the distribution is 0.043 s -1 , which gives a reasonable 
value of 0.25, when multiplied by A = 6.3 s. In fact, relax- 
ation of polymer molecules is a complex process, which 
involves a whole range of times. We measured relaxation 
of t p in our polymer solution after a sudden stop of a 
stationary shear flow with Wi = 6. The apparent re- 
laxation time, defined as (9Zn(r p ) / dt)^ 1 , increased from 
3 s to 40 s as the stress was decaying (see also Fig. 2b). 
Further, the PDF in Fig. 3 has pronounced exponential 
tails, so that flow events with large rate of extension oc- 
cur rather often. Therefore, we can conclude that the 
observed significant stretching of the polymers is quite 
consistent with the theoretical predictions. We believe 
that a similar kind of polymer stretching should occur 
in the usual, high Re, turbulent flows, when the velocity 
gradients are sufficiently large compared to A. 
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FIG. 1. Dependence of different parameters of flow on Wi. 
Wi was raised by 9%/min, 2.4%/min and 7%/min in the 
ranges [0.35, 2.1], [2.1, 8.7] and [8.7, 25], respectively, a) 
The ratio of the total stresses, cr/o~i am , curve 1. The ratio of 
the polymer contribution to the stress, t p /rf am , curve 2. b) 
Vrms in the center of the set-up, curve 1. Instrumental noise 
level in a laminar flow is shown for comparison, curve 2. c) 
Average first normal stress difference, Ni, in the turbulent 
flow, curve 1, and in a laminar flow between two plates, curve 
2. 

FIG. 2. Dependence of different parameters of flow on time, 
t. The time axis is compressed by a factor of 100 for t < 0. 
a) Flow velocity in the center corresponding to the curve 2 
below, b) The torque applied to the upper plate. Inset: curve 
2 in a vicinity of t = 0. 

FIG. 3. PDF for the longitudinal velocity gradients, 
dVj,/ (rd(j>) , estimated according to the Taylor hypothesis, 
with smoothing over about 1.3 mm. V<p was measured at 
z = 3.75 mm from the upper plate, at r = 15 mm, with time 
intervals At = 0.44 sec. The differences between the consec- 
utive velocity readings were divided by At and by V<p = 2.93 
mm/s. The thin line is a Gaussian fit. 
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